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Using a fermiology approach to the computation of the magnetic susceptibility measured by neu- 
tron scattering in hole-doped high-Tc superconductors, we estimate the effects on the incommensu- 
rate peaks caused by higher d-wave harmonics of the superconducting order parameter induced by 
underdoping. The input parameters for the Fermi surface and d-wave gap are taken directly from 
angle resolved photoemission (ARPES) experiments on Bi2Sr2CaCu20g+a; (Bi2212). We find that 
higher d-wave harmonics lower the momentum dependent spin gap at the incommensurate peaks 
as measured by the lowest spectral edge of the imaginary part in the frequency dependence of the 
magnetic susceptibility of Bi2212. This effect is robust whenever the fermiology approach captures 
the physics of high-T c superconductors. At energies above the resonance we observe diagonal incom- 
mensurate peaks. We show that the crossover from parallel incommensuration below the resonance 
energy to diagonal incommensuration above it is connected to the values and the degeneracies of 
the minima of the 2-particle energy continuum. 
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I. INTRODUCTION 

The imaginary part x"{oj,q) of the magnetic suscep- 
tibility x(a;,q) probed by neutron scattering in high- 
temperature (high-T c ) superconductors is characterized 
by a very rich dependence on energy and momentum 
transfer, temperature, and doping. 1-36 At temperatures 
well below the superconducting transition temperature 
T c and for a fixed frequency to (of the order of the d-wave 
gap maximum around optimal doping) within a finite 
range of frequencies, x" displays peaks at some symmetry 
related wave vectors in the Brillouin zone of the square 
lattice formed by planar Cu sites. The positions, heights, 
and widths of these peaks are temperature dependent 
and, in particular, are sensitive to the destruction of su- 
perconductivity above T c to a degree that depends on 
doping. As the energy transfered from the neutrons to 
the sample is varied, the positions in the Brillouin zone 
of the peaks in x" also vary (as well as their heights and 
widths), i.e., the peaks are dispersing. The detailed in- 
tensity distribution of the dynamical spin susceptibility 
depends on the high-T c superconducting family. 

For the YBa 2 Cu306+:r (YBCO) compounds it is ob- 
served that the separation in q space between four in- 
commensurate peaks at the symmetry related wave vec- 
tors q = (n ± 5, 7t) and q = (ir, tt ± S) decreases with 
increasing energy merging into a single resonance peak 
at the antiferromagnctic wave vector (tt, tt) and at an en- 
ergy oiAlmeV for optimally doped samples. 1-19 Cooling 
below the superconducting transition temperature opens 
up a doping-dependent spin gap, which is proportional 
toT c . 15 

In the La2- 2; Sr 2 ;Cu04 (LSCO) family the dynamical 
magnetic susceptibility follows a similar dispersion. At 



energies above a spin gap four incommensurate peaks on 
the horizontal and vertical lines passing through (tt, tt) 
in the Brillouin zone appear and disperse towards the 
antifcrromagnetic zone center as the energy transfer is 
increased. 20-31 Although the incommensurate peaks in 
LSCO seem to join at the commensurate wave vec- 
tor q = (tt,tt), no resonance peak has been observed in 
this compound to this date. Another remarkable obser- 
vation in LSCO is that a spin gap has not been observed 
in the extreme underdoped and overdoped regimes. 29 ' 30 

On Bi 2 Sr 2 CaCu 2 8+2; (Bi2212) and Tl 2 Ba 2 Cu0 6+:l; 
(TBCO) samples only a few inelastic neutron scatter- 
ing measurements have been performed so far. For both 
compounds a resonance peak was observed, 32-36 but an 
incommensurate signal below the resonance has not yet 
been measured, perhaps due to the limited size of single- 
crystals presently available. 

In this paper we will take the point of view that incom- 
mensurate and commensurate peaks of x" have a com- 
mon origin as they appear to be smoothly connected in 
the YBCO family. 9 ' 11 ' 12 ' 15 ' 37 This hypothesis seems hard 
to reconciliate with theoretical scenarios based on the 
existence of dynamical stripes or on the SO (5) approach 
to high-T c superconductivity. In the dynamical stripe 
scenario incommensurate peaks are the natural descen- 
dants of the static charge and spin long-range order seen 
in La 2 _ a; Ba :i; Cu04 at x = 1/8, 38 say. 39-45 In the SO(5) 
scenario a resonance at (tt, tt) appears naturally as a re- 
sult of a bound state in the spin-triplet particle-particle 
channel. 46-52 Commensurate and incommensurate peaks 
of x" are smoothly connected in a scenario in which it is 
assumed that strongly renormalized quasiparticles close 
to the Fermi surface interact with a residual on-site Hub- 
bard or nearest-neighbor antifcrromagnetic interaction, 
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in short a fcrmiology approach. However it remains to 
be shown if any disagreement between the line shapes 
for the magnetic susceptibility predicted by fcrmiology 
and the measured line shapes is a mere quantitative one 
that results from the inherent approximations involved in 
fcrmiology or is of a more qualitative nature that reflects 
physics beyond the reach of fcrmiology. 

Common to all fcrmiology scenarios 53-102 (see also 
Refs. 103-105 for related works on the "bare" magnetic 
susceptibility) is the random phase approximation (RPA) 

( \ q) n r , 

XHPA(-,q)= 1 + Mwjq)nKq) (1-1) 

to the magnetic susceptibility. Here, n(w, q) is 
the ("bare") magnetic susceptibility of noninteracting 
fcrmionic (BCS) quasiparticles that depends sensitively 
on the Fermi surface (the superconducting BCS gap) 
above (below) T c . The functions g(ui,q) and h(oj,q) 
are model dependent: A RPA treatment of the single- 
band Hubbard model with on-site repulsion U yields 
g{u,q) = n(w,q) and h(u,q) = -?7. 53-63 A RPA 
treatment of a single-band of fermionic (BCS) quasi- 
particles with a residual interaction such as a nearest- 
neighbor antiferromagnetic interaction of strength pro- 
portional to J yields g(u>,q) = H(uj,q) and h(oj,q) = 
(J/2)(cosq x + cos q y ), 64-67 as is often done in the slave- 
boson treatment 68-80 or the l/z expansion with z the 
number of nearest-neighbors of the t — J model. 81-84 A 
RPA treatment on the bare propagator Xo( w ;Q) of col- 
lective spin-1 excitations due to a weak coupling as mea- 
sured by the coupling constant g with otherwise nonin- 
teracting fermionic (BCS) quasiparticles yields g(ui, q) = 
X (w,q) and h(uj,q) = -g 2 x (w, q). 85-89 

The RPA approximation (1.1) has been improved in 
three ways. First, the single-band Hubbard model can 
be generalized to the three-band Hubbard model. 90-94 
Second, the feedback effect of the magnetic fluctua- 
tions encoded by Eq. (1.1) on the propagator of the 
fcrmionic (BCS) quasiparticles can be included self- 
consistcntly through the so-called fluctuation exchange 
(FLEX) approximation. 95-100 Third, all unaccounted for 
interactions among the fermionic (BCS) quasiparticles 
can be included in a phcnomcnological way by substi- 
tuting in Eq. (1.1) u) by oj + iT with T a positive function 
of uj, temperature T, doping x, and high-T c supercon- 
ducting family. 101 ' 102 

The goal of this paper is to assess the effects on 
the incommensurate peaks of the magnetic susceptibil- 
ity x caused by the presence of a higher <i-wave har- 
monic in the superconducting order parameter of un- 
derdoped Bi2212. Indeed, it is observed in several an- 
gle resolved photoemission (ARPES) experiments on un- 
derdoped Bi2212 that there are deviations from a pure 
d-wave order parameter. 106-108 These deviations result 
in a rounding of the superconducting gap in the vicinity 
of its nodes. 109 We interpret these deviations as the sig- 
nature of a d-wave order parameter extending to next- 
nearest- neighbor bonds between planar Cu sites, i.e., 



they push closer to the anisotropy ratio v A /v F between 
the slope of the superconducting gap v A and the Fermi 
velocity v F , which are tangent and perpendicular to the 
Fermi surface, respectively. 

Below, we adopt a phenomenological fermiology sce- 
nario by which wc compute the RPA susceptibility (1.1) 
assuming that IL(uj, q) is calculated with the BCS dis- 
persion measured by ARPES while g(oj, q) = n(w, q) 
and h(ui, q) = —U. To simplify the matter and to isolate 
the effect of decreasing v A {x)/v F (x) with underdoping, 
i.e., decreasing x, we treat residual lifetime effects on the 
BCS quasiparticles with the substitution 

uj — >uj + iT. (1.2) 

For reference T = 1 meV is of the order of the energy res- 
olution of most experiments. The lowest spectral edge in 
Xrpa( w ' l) considered as a function of uj with q held fixed 
defines the q-dependent spin gap. We find that increas- 
ing the higher harmonics by an increment of 15% relative 
to optimal doping, which corresponds to a decrease of T c 
by 28% in Bi2212 as was observed in Ref. 107, decreases 
the q-dependent spin gap by an amount that depends on 
the wave vector (see Figs. 3 and 4) and agrees qualita- 
tively with YBCO measurements. This result is a robust 
feature of all fcrmiology scenarios and is consistent with 
the experimental observation of a fast decreasing (or even 
vanishing) spin gap with underdoping. 15,29 

The paper is organized as follows. The fcrmiology ap- 
proach is defined through a RPA magnetic susceptibility 
in Sec. II. Numerical results are presented in Sec. Ill 
and interpreted in Sec. IV. Wc close with a summary in 
Sec. V. 



II. RPA MAGNETIC SUSCEPTIBILITY 

Perhaps the most important lesson that can be inferred 
from ARPES data is that the low energy excitations in 
the superconducting state of Bi2212 arc sharply defined 
quasiparticles obeying the BCS dispersion 110 

= V^k + A l (2-la) 

in the close vicinity to the Fermi surface defined by the 
condition 

£k = (2.1b) 

and with a superconducting gap consistent with a <i-wave 
symmetry, i.e., 

A k = — A k , , (k' x7 k y ) = (±k y ,^pk x ) , 

(2.1c) 

A k = A k „ (k' x ,k' y ) = (±k x ,Tk y ) . 

This fact suggests that it might be plausible to treat the 
superconducting state of high-T c superconductors as a 
conventional BCS superconductor in the close vicinity 
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to the Fermi surface, an assumption that we will make 
from now on. Following Norman in Ref. 62 the measured 
Fermi surface is fitted from the tight-binding expansion 
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(2.1d) 



with a n 



and a , b a pair of orthogonal vectors 



joining a site of the square lattice to two of its j-th nearest 
neighbors. Similarly, the d-wave BCS gap is fitted from 
the tight-binding expansion 



A 1 (cos k x — cos k y ) 

- cos 2k y 



+A 2 (cos 2^ — cos2fcj,) (2-le) 

+A3 (cos 2k x cos k y — cos k x cos 2k y ) 
_i_ 



In practice we truncate the expansions at j = 5 for the 
Fermi surface and to the first two terms for the d-wave 
BCS gap. As opposed to the slave-boson approach which 
attempts to capture the doping dependence of the i^'s 
and the A^-'s from the large ?7-limit of the Hubbard 
model, we will take the phenomenological point of view 
that the BCS dispersion (2.1) of the "bare" quasiparti- 
cles is an input deduced from ARPES data for the Bi2212 
family. 

We will assume in this paper that anomalous features 
at energy scales far away from the Fermi energy, say 
of the order 2max k A k such as is the case for the in- 
commensurate and commensurate peaks observed at op- 
timal doping with inelastic neutron scattering, can be ac- 
counted for by postulating the existence of "bare" nonin- 
teracting BCS quasiparticlcs whose dispersion is given by 
a fit to ARPES data and which interact weakly through 
some small residual interactions which we take to be an 
on-site Hubbard repulsion with coupling constant U . In 
this spirit, the renormalization of the "bare" magnetic 
susceptibility 
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FIG. 1: Plot of the real parts XkpA(0,q) (a) and ll'(0,q) (b) 
as a function of q in the upper-right quadrant of the Brillouin 
zone for the BCS dispersion of Bi2212 at optimal doping (see 
table I), (c) Plot of II" (cj, tt, n)/[nv 2 (uj, n, tt)] as a function of 
lj > in the relevant energy range for the BCS dispersion of 
Bi2212 at optimal doping (see table I). The inset displays the 
uj dependence of the imaginary part II" (w, n, n) (thick line) 
and ni>2(u}, n, tt) (dashed line). The temperature is here taken 
to be T = OK while a damping F — O.lmeV is used. 
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is, within the RPA, 

XRPA(w,q) 
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l-UIL((j,q)- 



(2.2b) 



Here, is the number of sites on a square lattice with 
periodic boundary conditions, f(x) = [exp(/3x) + 1] _1 is 
the Fermi distribution (/3 the inverse temperature), and 
the BCS dispersion (2.1) is used. 

In Fig. la and Fig. lb we present a plot of the static 
limit of the RPA magnetic susceptibility (2.2b) and of the 



static "bare" magnetic susceptibility (2.2a), respectively. 
It is seen that there are no significant differences between 
the static limits of the "bare" and RPA magnetic suscep- 
tibilities for small momentum transfer q (small forward 
scattering) but that for large momentum transfer, q close 
to (tt, 7r), the renormalization effects are of order 1 as they 
signal the proximity to an antifcrromagnctic instability 
of the BCS ground state. These renormalization effects 
modify the quasiparticle self-energy and have been pro- 
posed as an explanation to the peak-dip-hump shape seen 
in the ARPES signal as a function of frequency. 74 ' 112 In 
this paper we do not consider this feedback on the quasi- 
particle self-energy. However, since the effects of higher 
d-w&ve harmonics on the RPA magnetic susceptibility 
are robust, we expect them to survive in a selfconsistent 
approach. 



We close this section by noting that the imaginary part 
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TABLE I: Values of the hopping and gap parameters used to 
fit the Fermi surface (as in Ref. 62) and gap (see Ref. 107) 
in the Bi2212 family measured with ARPES as a function 
of doping. The gap parameters for T c = 87K, T c = 837^ 
and T c = 75K follow from a fit to ARPES data (Ref. 107), 
and for T c = 68K from an extrapolation. The doping de- 
creases when reading the columns from left to right: The 
second, third, fourth, and fifth columns correspond to opti- 
mally doped, slightly underdoped, and underdoped samples, 
respectively, as indicated by their decreasing T c . The critical 
temperatures are in units of Kelvins while all other quantities 
are in units of meV. 

Parameters T c = 87K T c = 83K T c = 75K T c = 68K 
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U"(ui,q) of the "bare" susceptibility (2.2a) simplifies to 

n >«i) ^E c qV%- £ 2(q,k)), 

k 

^ir^(i- ' t+ 'Xl q A r Ak )' (23) 

£ 2 (q,k) =£ k+q + £ k , 

in the zero temperature limit and for positive frequencies 
as is relevant in inelastic neutron scattering. Wc shall sec 
in Sec. IV that the numerical results from Sec. Ill can be 
understood in terms of the frequency dependence of the 
2-particle density of states (DOS) 

"a(^q):=2^E*( w -- B 2(q.k)) ( 2 -4) 

k 

since it differs very little from that of n"(w,q)/7r for a 
fixed wave vector q close to the commensurate vector 
(■7T, 7r) and at very low temperatures as is illustrated in 
Fig. lc. 

III. NUMERICAL RESULTS 

We have calculated numerically the RPA (2.2b) to the 
magnetic susceptibility with the inclusion of lifetime ef- 
fects implemented through substitution (1.2). Hereby, we 
are using the BCS dispersion (2.1) with the hopping and 
superconducting parameters extracted from fits to the 



Fermi surface and gap measured with ARPES within the 
Bi2212 family for the underdoped and optimally doped 
regimes (see Refs. 62,107,113-115). These parameters 
can be found in table I. For each doping concentration 
the value of the on-site Hubbard repulsion U is chosen so 
that the energy of the resonance is consistent with unpo- 
larized neutron data in Refs. 32,34 and with the indirect 
determination of the resonance energy via the peak-dip- 
hump feature from APRES data in Refs. 116,117. The 
value of the coupling constant U increases with under- 
doping, a trend consistent with the naive expectation 
that moving away from half-filling in the Hubbard model 
reduces the effect of the strong interactions. ARPES ex- 
periments on Bi2212 indicate that the Fermi surface is (if 
at all) only weakly dependent on doping. 114,115 To sim- 
plify matters, wc have implemented this observation by 
keeping the hopping parameters unchanged as a function 
of doping. 

To perform the summations in formula (2.2b), we 
meshed the Brillouin zone with 1024 x 1024 points. In 
order to reduce the effects of a finite lifetime in Fig. 1, we 
reduced T in Eq. (1.2) from the experimental resolution 
~ lmeV to O.lmeF. The temperature is taken to be 
vanishing in most instances except in Fig. 3 and Fig. 6 
where T = 5K. 

It is convenient to use polar coordinates to represent 
the wave-vector dependence of the superconducting gap 
on the Fermi surface. To this end, define along the arc of 
the Fermi surface (2.1b) that belongs to the upper- right 
quadrant of the Brillouin zone the angle < < 7r/2 
through (see Fig. 2a) 

<j) = arctan ( - - ky ) . (3.1) 

V 7T - k x J 

In terms of this angle, define the auxiliary gap function 

A(0) := A [B cos(20) + (1 - B) cos(60)] , 

A := max A k . (3.2) 

k^Formi surface 

The auxiliary gap function A(0) obeys the same ci-wave 
symmetry as the gap function A k does on the Fermi 
surface. The parameter B entering A(<^) is taken to 
be a function of doping x that can be determined for 
Bi2212 from ARPES measurements (see table I). 107 The 
agreement between A(<^) and A k for k on the Fermi sur- 
face is excellent. The Fermi velocity v A tangent to the 
Fermi surface at the node of the gap is proportional to 
(dA) j (d<f>) at (j> = it/ A. It decreases with undcrdoping as 
is apparent from table I. The gap function A k on the 
Fermi surface is depicted in Fig. 2c for the parameters of 
table I. 

In Figs. 3a-c we illustrate for three different wave vec- 
tors depicted in Fig. 2b, q = (ir, tt), the antiferromag- 
netic wave vector, q = (1.267T, 7r), a wave vector on the 
horizontal symmetry axis passing through (tt, it), and 
q = (1.065-7T, 0.8057r), a wave vector off the symmetry 
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FIG. 2: Panel (a) displays the Fermi surface of Bi2212 with 
to, . . . ,t& as in table I. The angle cj> measures the position 
along the Fermi arc centered about (tt, n) (hole-like) in the 
Brillouin zone. Panel (b) displays the reciprocal q space to- 
gether with the symmetry axes (parallel and diagonal) and 
the three wave vectors (n, ir), (1.267T, 7r), and (1.0657T, 0.8057r). 
Panel (c) displays the absolute value of the gap function 
| A (0)| on the Fermi surface as a function of Fermi surface 
angle cj> for Bi2212 with the parameters of table I. Plot- 
ted are the gaps for optimally doped samples (solid line), 
slightly underdoped (long dashed line) , and underdoped sam- 
ples (dashed and dotted lines). For the momentum transfer 
q = (1.26-7T, 7r) vertical lines at the Fermi surface angles, where 
the pairwise nested wave vectors ki , ki + q, k2 and k2 + q lie, 
are drawn (see Sec. IV). 



axes, how the dependences on frequency of the imagi- 
nary part Xrpa + £T, q) of the RPA magnetic suscepti- 
bility (2.2b) changes with undcrdoping. In each panel 
the four curves correspond to the four dopings in ta- 
ble I. As the doping decreases (i.e., decreasing T c ) the 
leading edge moves to the left (i.e., to lower frequen- 
cies). The q-dcpendcnt spin gap, which we define by 
the lowest spectral edge in Xrpa considered as a func- 
tion of frequency, decreases with decreasing doping for 
both the commensurate wave vector (tt, 7t) and the in- 
commensurate points (1.26tt, tt) and (1.065tt, 0.8057r) on 
the symmetry axis and off the symmetry axis passing 
through (tt, tt), respectively. This behavior is consistent 
with the observed doping dependence of the q-depcndent 
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FIG. 3: Panel (a) displays the frequency dependences of 
Xrpa^iQ) a t the commensurate wave vector q = (n,n) for 
four different values x of doping in the Bi2212 family with 
the BCS dispersions taken from table I. Curves with the 
leading edge moving to the left (i.e., to lower frequencies) 
have decreasing doping x (T c ). Panels (b) and (c) are the 
same as panel (a) except for the incommensurate wave vectors 
q = (1.26tt, tt) and q = (1.065tt, 0.805tt) being held fixed. For 
all three wave vectors the q-dependent spin gap decreases with 
underdoping. The temperature is here taken to be T = 5K 
while a damping V = lmeV is used. 



spin gap in the LSCO and YBCO families. The intensi- 
ties of the peaks in all panels of Fig. 3a-c increase with 
underdoping. This increase in intensity is accompanied 
by a narrowing of the line shape in panels b and c. 

In Fig. 4 the intensities of the imaginary part of the 
RPA spin susceptibility are shown as a function of frc- 
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quency and wave vector for the four doping concentra- 
tions of table I. At the corresponding resonance ener- 
gies (39meV, 37meV, 34meV, and 31meV, respectively) 
the peaks are at the (tt, tt) point. When the frequency 
is reduced from the resonance energy, dominant incom- 
mensurate peaks in the (tt + S, tt) and (tt, tt + d) direc- 
tion occur, and a subdominant structure on the diagonal 
lines which pass through two incommensurate points in 
q space shows up. The peaks arc dispersing with a down- 
ward curvature. Above the resonance energy the domi- 
nant peaks cross over from the parallel to the diagonal 
symmetry axes passing through (tt,tt). The same hap- 
pens below 19.2meV (T c = 87K), 17.6meV (T c = 83K), 
15.2meV (T c = 75K), and 12.8meV (T c = 68K), re- 
spectively. (We will then speak of parallel and diago- 
nal peaks, respectively.) For frequencies larger than the 
resonance energy incommensurate peaks on the diago- 
nal symmetry axes passing through the antiferromagnetic 
wave vector (tt, tt) have recently been observed in both 
La2-a;Ba K Cu04 38 and YBCO. 19 The crossover to dom- 
inant diagonal peaks at low energies has not been ob- 
served in experiments. If the wave vector is restricted to 
the horizontal symmetry axes in the Brillouin zone, the 
spin gap A SG (q x , tt) remains finite for all values of q x . 
However, if the wave vector lies on the diagonal symme- 
try axes the spin gap A SG (q x , q x ) vanishes for certain 
values of q x = q . 

The dispersion of the peak positions of Xrpa f° r the 
four doping concentrations of table I and with q moving 
along the horizontal or diagonal lines passing through 
(tt, tt) arc shown in Fig. 5a and Fig. 5b, respectively. 
With undcrdoping the bell-like shape of the peak disper- 
sion moves to lower energies implying a decrease of the 
q-dependent spin gap (see also Figs. 3 and 4). The de- 
gree of incommensuration measured by the separation 26 
of the two peaks of Xrpa decreases as one moves from the 
optimally doped to the undcrdoped regimes. This agrees 
with neutron scattering measurements on the YBCO and 
LSCO families. 15 ' 27 



IV. DISCUSSION 
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FIG. 4: Panels (a)-(d) display the intensities of Xrpa^i fe, tt) 
in the (uj,q x ) plane at fixed q y — n for four different doping 
concentrations in the Bi2212 family with the BCS dispersions 
taken from table I. Panels (a), (b), (c), and (d) correspond to 
samples with a transition temperature T c of 87K, 83K, 75K, 
and 68^, respectively. As a guide to the eye contour lines are 
drawn at an intensity which is about 0.5% of the maximum 
intensity of the dynamical spin susceptibility. Panels (e)-(h) 
are the same as panels (a)-(d) except for the q-vector which 
runs along the diagonal instead of the horizontal line in re- 
ciprocal space. The temperature is here taken to be T — OK 
while a damping T = lmeV is used. 



In this section we review the mechanism causing com- 
mensurate resonance peaks and incommensurate peaks 
in the numerical simulations of the imaginary part of the 
RPA magnetic susceptibility. The discussion is kept brief 
as it has been thoroughly described in the fermiology pa- 
pers of the early 90's with the exception of one aspect 
based on some point group symmetries of the BCS dis- 
persion. 

Up to multiplication with the Bosc-Einstcin distribu- 
tion, inelastic neutron scattering has direct access to the 
imaginary part x"(o^c[) of the magnetic susceptibility 



x(w,q), which is given by 

/, , v n"(^, q ) 

[l-[/nV,q)] 2 +[im"(^,q)] 2 

(4.1) 

within the RPA approximation [x'(w, q) and II'(w,q) 
denote the real parts of x(^,q) and U(ui,q), respec- 
tively]. The "bare" imaginary part II" (w, q) of the 
magnetic susceptibility controls x^ pA (w,q) in two ways. 
First, Xrpa( w 'CI) vanishes whenever n"(w,q) does and 
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1 ^ UH'(u>,q). Second, at q held fixed, any steplikc 
discontinuity in the frequency dependence of II" (w, q) at 
some frequency w„(q) results in a logarithmic divergence 
of the frequency dependence of II' (w, q) at w n (q) through 
the Kramers-Kronig relation between the imaginary and 
real parts of causal response functions. This in turn guar- 
antees that (i) the dynamical Stoner condition 



1 - UU'(u},q) = 



can be met at the frequency 



<(q) < w «(q) 



(4.2a) 



(4.2b) 



and (ii) x^ PA (oj,q) acquires a pole at wj(q) and peaks 
at o>*(q), where n > 1 indexes the remaining jump dis- 
continuities of II" (w, q). 

For any small but finite damping T, II" (co + iT, q) be- 
comes a nonvanishing and continuous function of u> > 
when holding q fixed. Similarly, the logarithmic singu- 
larities in the cu dependences of II' (w, q) are cutoff by a 
finite r under the substitution (1.2). Under this substi- 
tution the dynamical Stoner criterion (4.2a) can only be 
met for a sufficiently large size of the step in II" (w, q), 
and the pole at ojI in Xp PA (u;,q) turns into a peak of 
finite height 



XrpaK +«r,q) 



U 2 W{loI +«r,q) 

with the full width at half maximum (FWHM) 



FWHM 



2n"(w + ir,q) 



(dU'(uj + iT, q)/du) 



(4.2c) 



(4.2d) 



> 

CD 

E 



40 



30 



20 



10 




T T C =87K 

O T c =83 K 

♦ T C =75K 

V T =68 K 




a) 



0.6 



0.8 



1.0 

q x m 



1.2 



1.4 



40 



^ 30 
> 

CD 



3 20 



10- 



b) 




▼ T =87K 

O T G =83 K 

♦ T =75K 

V T =68 K 



1 1 1 

0.7 0.8 0.9 1.0 1.1 

(q.q) (i) 




1 .2 



1.3 



To sum up, the condition for dynamical Stoner en- 
hancement of Xp PA (tj,q), a collective effect, can be re- 
duced to the condition for II" (u, q) to have steplikc dis- 
continuities as a function of frequency when holding the 
wave vector fixed. 

In view of Eqs. (2.3) and (2.4) steplike discontinuities 
in the frequency dependence of II" (w, q) with q held fixed 
(at zero temperature and with the inverse lifetime T in- 
finitesimally small) are closely related to steplikc discon- 
tinuities in the frequency dependence of the 2-particle 
DOS (2.4) at fixed q. These jump discontinuities in 
v 2 (io, q) occur whenever the 2-particle energy E 2 (q, k) in 
Eq. (2.3) reaches a local minimum cu n at the q-dependent 
wave vector k ?l , in which case the size of the step is pro- 
portional to 118 



1/2 



(4.3) 



< 



d 2 E 2 (q,k) 



m 



(L 2 ) dk 12 dk u 



The steplikc discontinuity in the 2-particle DOS z/ 2 (a;,q) 
at to n is present in II" (a;, q) if and only if the coherence 



FIG. 5: (a) Peak positions of Xrpa( lj i c i) f° r fixed q y — n as 
a function of q x and energy with the BCS dispersions taken 
from table I for four different samples with a transition tem- 
perature T c of 87 K (solid triangles), 83K (open diamonds), 
75 K (solid diamonds), and 68 A" (open triangles), (b) Peak 
positions of Xrpa( w i q) as function of frequency and momen- 
tum transfer with q along the diagonal q x = q y for the same 
four doping concentrations as in panel (a). Above the reso- 
nance the diagonal peaks are the dominant structure. The 
same happens at very low energies. At intermediate energies, 
the diagonal peaks are subdominant with the parallel peaks 
being dominant (see Fig. 4). The temperature is here taken 
to be T — OK while a damping F — lmeV is used. 



factor ^ from Eq. (2.3) is non-vanishing in a neigh- 
borhood of k n . The steplikc jump is turned into a smooth 
increase otherwise. 104 

It is shown in appendix A that, when q is chosen not to 
be on any of the diagonals of the Brillouin zone or the hor- 
izontal and vertical lines passing through the commensu- 
rate point (-7T, 7r) (see Fig. 2b), there are four distinct but 
2-fold degenerate local minima uj 1 < uj 2 < lu 3 < uj^ of 
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E 2 (q, k) for the parameters of table I located at the q- 
dependent wave vector k 1; k 2 , k 3 , and k.4, respectively, 
in the Brillouin zone. When q is chosen to be on the diag- 
onals that pass through the commensurate point (tt,tt), 
but away from it, there are are three distinct local min- 
ima u>i < uj 2 < uj 3 at ki, k2, and k3, respectively, which 
are either doubly degenerate (lo 1 and u> 3 ) or 4-fold degen- 
erate (w 2 ). If the momentum transfer lies on the horizon- 
tal or vertical lines passing through the commensurate 
point, but away from it, there are two distinct but 4- fold 
degenerate local minima uj 1 < uj 2 of E 2 (q, k) located at 
the q-dependent wave vector ki and k 2 , respectively. Fi- 
nally, all local minima collapse to one global but 8-fold 
degenerate minimum lo 1 located at k x when q = (tt, tt) 
or one of its symmetry related images. 

Two criteria control the size of the first step with in- 
creasing lo: (i) How large the effective masses rnf^' 2 ' are, 
i.e., how flat the 2-particle dispersion is, (ii) how large 
the degeneracy of the global minimum is. It turns out 
that, when q is chosen to be at the commensurate anti- 
ferromagnetic vector (tt, tt), the geometrical mean of the 

effective masses m^ 1 ' 2 ' becomes very large because of the 
proximity of k x to an extended saddle-point for the pa- 
rameters of table I. Furthermore this effect is magni- 
fied by the 8-fold degeneracy of the global minimum of 
the 2-particle dispersion at q = (tt, tt). Taken together 
conditions (i) and (ii) insure that the dynamical Stoncr 
criterion is met at an energy 

< = ( 4 - 4 ) 

sufficiently far below the 2-particle continuum threshold 
for the resonant nature of the resulting spin-1 collective 
excitation not to be washed out by a finite temperature 
of 5 if or a damping of 1 meV. The 8-fold degeneracy 
is partially broken to a 4-fold degeneracy when q moves 
away from (tt, tt) along the horizontal or vertical lines 
passing through (tt, tt) (parallel incommensuration) with 
kj moving towards the nodal line and k 2 moving away 
from it (see Fig. 2c). For energies below the commensu- 
rate resonance energy 

w * = w pc < w af> ( 4 -5) 

this turns the resonant commensurate peak into a down- 
ward dispersing incommensurate peak, which eventu- 
ally becomes nonresonant due to the failure to meet 
the Stoner criterion. By energy conservation the dom- 
inant incommensurate peaks cross over to the diagonals 
of the Brillouin zone at frequencies much lower than lu^ f 
since q then connects regions close to the nodes of the 
BCS dispersion. 103 When the wave vector is pushed away 
from (tt,tt) along the diagonals, the 8-fold global mini- 
mum splits into three distinct minima. The first min- 
imum lui being lowest in energy leads to the diagonal 
peaks at low energies. Together with the lowest mini- 
mum of the 2-particlc energy for q off the symmetry axes 
it is also the cause of the subdominant structure on the 
diagonal lines passing through (tt + S, tt) and (tt, tt + S), 



say, at intermediate energies. The second 4-fold degen- 
erate minimum lies in general at energies just above the 
resonance w^ F , but below the 8-fold degenerate minima 
at q = (tt, tt), and is thus responsible for the crossover 
to diagonal incommensurate peaks above the resonance 
energy. Finally, by choosing q away from the symme- 
try axes one spreads the spectral weight of the commen- 
surate resonance among four distinct 2-fold degenerate 
local minima, three of which are in the 2-particle contin- 
uum thereby loosing the most in height and sharpness of 
the peak relative to the line shape of the commensurate 
resonance as a function of frequency. 

To conclude, the crossover from diagonal to parallel 
incommensuration with increasing frequency at low en- 
ergies and from parallel to diagonal peaks with increasing 
frequency above the resonance energy is brought about 
by the onset of 4-fold degenerate minima of the 2-particle 
dispersion E^q, k) at the corresponding energies. This 
effect is quite sensitive to both the shape of the BCS dis- 
persions and the details of the residual interaction among 
the BCS quasiparticles: On the one hand, an antifcrro- 
magnetic interaction with h(ui, q) = (J/2)(cosq x +cosq y ) 
leads to an enhancement of the diagonal peaks com- 
pared to the parallel ones, which potentially destroys the 
crossover to parallel incommensuration upon increasing 
frequency towards the resonance. 62 On the other hand, 
it was observed in Refs. 76,78 that the parallel peaks re- 
main dominant at energies below but not too far from 
the resonance in a slave-boson mean-field approach to 
the t — t' — J model. 

At low energies the intensity of Xrpa( w + q) for a 
fixed wave vector dies off quickly (see Figs. 3 and 4). In 
the following we shall define the lowest, sharp spectral 
edge in x'kpa^^ l) considered as a function of frequency 
at a fixed wave vector q to be the q-dependent spin gap. 
By Eqs. (2.3) and (4.1), the q-dependent spin gap for 
vanishing damping T is given by the threshold to the 
2-particle continuum provided there is no resonant exci- 
tation below the threshold. Otherwise, the q-dependent 
spin gap is determined by the energy of the lowest bound 
state. Any sharp spectral edge in Xrpa^'I) turns into 
a smooth edge by the inclusion of finite lifetime effects 
through substitution (1.2) in which case, to a first ap- 
proximation, the q-dependent spin gap is given by the 
position of the half maximum of the leading edge on the 
low frequency side. 

The effects on the q-dependent spin gap caused by the 
change of the gap parameters in table I induced by under- 
doping can be best understood by use of arguments based 
on pseudo Fermi nesting. A good approximation for the 
position of the q-dependent (local) minima of E 2 (q, k) is 
obtained by requesting that both k n and k n + q lie on 
the Fermi surface 

= £ k „ = e k „+ q ( 4 -6a) 
and that the coherence factor ^~ be maximal, i.e., 

sgn(A kn A kn+q ) =-1. (4.6b) 
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For a fixed wave vector q the steplike discontinuity of the 
frequency dependence of II" (to, q) is then approximately 
located at the energy 



A(kJ + A(k n + q) 



(4.6c) 



where k n are the solutions of the equation set (4.6a). In 
particular, uj 1 determines approximately the threshold to 
the 2-particle continuum. 

At frequencies much smaller than the resonance en- 
ergy the approximate 2-particle threshold ui\ is ef- 
fectively governed by excitations whose pscudo-Fcrmi- 
nesting vectors connect parts of the Fermi surface which 
are close to the nodes. Hence, at these energies the de- 
creasing slope of the gap function (dA)/(d(j)) at <f) = 7r/4 
results in a decrease of the approximate 2-particlc thresh- 
old uji (see vertical lines ki and ki + q in Fig. 2c). 

For frequencies which are comparable to the resonance 
energy however, the 2-particle threshold increases with 
the change of the gap parameters induced by undcrdop- 
ing. This can be understood from the fact that at fre- 
quencies close to, or higher than the resonance energy 
the approximate 2-particle threshold uj\ is controlled by 
excitations with pseudo-Fcrmi-ncsting vectors that con- 
nect parts of the Fermi surface which are in the vicin- 
ity of (0, 7r) or one of its symmetry related images. In 
these parts of the Fermi surface the increase of the gap 
maximum Ao causes the gap function to increase with 
underdoping (see vertical lines k2 and k2 + q in Fig. 2c) . 

On top of the effects caused by the change of the 
gap parameters there is an additional modification 
of XrpaO^ + q) with underdoping induced by the in- 
crease of the coupling constant U. Owing to the failure 
to meet the Stoner criterion, the frequency of the in- 
commensurate peaks at a given wave vector q and at 
energies much smaller than the resonance energy is not 
changed by a not-too-large increase in U. Thus, in the 
limit of a T — 0, the q-dependent spin gap is unaffected 
in this energy range by a not-too-large increase of the 
coupling constant. For energies close to the resonance 
energy though, the increase of U with underdoping 
causes both the q-dependent spin gap and the frequency 
of the (in)commensurate peaks to decrease. The quanti- 
tative role played by U can be illustrated by reproducing 
the frequency dependences of the RPA susceptibility in 
Fig. 3 with all BCS parameters unchanged but with U 
fixed to its value at optimal doping. The result is dis- 
played in Fig. 6a and Fig. 6b where one sees that: (i) 
The resonance energy moves with underdoping to higher 
energies as a result of the increasing gap maximum in 
panel (a), and (ii) the leading edge moves with under- 
doping to lower energies for the wave vector (1.267T, ir) in 
panel (b). (The position of the leading edge is defined 
to be at the half maximum of the edge on the low fre- 
quency side.) We see that the spin gap value at a wave 
vector sufficiently far away from (ir, tt) is rather insensi- 
tive to keeping U constant as in Fig. 6b or to adjusting 
U so as to get the correct resonance energy as in Fig. 3b. 




10 20 30 40 50 60 
co [meV] 



FIG. 6: Panel (a) displays the frequency dependences of 
Xrpa^iQ) a t the commensurate wave vector q = (n,n) for 
four different values x of doping in the Bi2212 family with 
the BCS dispersions taken from table I but with U fixed 
to 165meV for all doping concentrations. Panel (b) is the 
same as panel (a) except for the incommensurate wave vec- 
tor q = (1.267r,7r) being held fixed. Panel (c) displays the 
frequency dependences of XRPA^q) a t the incommensurate 
wave vector q = (1.267T, 7r) for four different values x of doping 
in the Bi2212 family without inclusion of higher harmonics, 
i.e., with vanishing A 2 . As explained in the text we take here 
Ai to be 18.8meV, !9.9meV, 20.5meV r , and 22.1meV , for 
samples with a T c of 87K , 83K, 75K, and 68K, respectively, 
and U increases with underdoping from 170meV at optimal 
doping to 184meV , 196meV, and 21ImeV in the underdoped 
regime. In this way, U is adjusted to reproduce the positions 
of the resonance as observed in Bi22I2. The temperature is 
here taken to be T — 5K while a damping T = lmeV is used. 
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We conclude that the dominant doping dependence of the 
q-dependent spin gap comes from the doping dependence 
of the higher harmonics for a wave vector sufficiently far 
away from (tt, tt). 

This insensitivity depends in a crucial way on how large 
U becomes with underdoping. This can be illustrated in 
a rather dramatic way by switching off the BCS parame- 
ter A 2 for all dopings. The parameter A x is then chosen 
so that the maximum of the superconducting gap on the 
Fermi surface A agrees with the values in Ref. 107, i.e., 
we find that A 1 equals 18.8meV, 19.9meV, 2Q.5meV, 
and 22.lmeV, for samples with a T c of 87 K, 83K, 75K, 
and 68K, respectively. For the coupling constant U we 
demand that its value is chosen so as to reproduce the 
same resonances as in Fig. 3a, i.e., we find that U takes 
the values 170meV, 184meV, 196meV, and 211meV, 
when T c takes the values 87K, 83K, 75K, and 68K, re- 
spectively. Evidently, the increase of U with underdoping 
is now much stronger than in table I. Having recalibrated 
U to the new parameters of the BCS dispersion, we plot 
in Fig. 6c the frequency dependence of x"(w + iT, q) for 
the wave vector (1.267T, 7r). Comparison of Fig. 6c with 
Fig. 3b shows that a broad peak at optimal doping can 
be turned into a resonance with underdoping due to a 
too strong increase in U. 

To summarize, we expect that the q-dependent spin 
gap always decreases with underdoping. Far away from 
the antifcrromagnetic wave vector (tt, tt) this is mostly a 
consequence of the decreasing slope of the gap function at 
the node, whereas close to (it, tt), it is a result of both the 
increasing coupling constant U and the increasing higher 
harmonics. This is confirmed by the numerical simula- 
tions presented in Sec. III. The effect of a decreasing 
q-dependent spin gap upon changing the gap parame- 
ters induced by underdoping is present in all fermiology 
scenarios. For example, we have also performed calcula- 
tions with a residual nearest-neighbor antifcrromagnetic 
interaction instead of the on-site Hubbard repulsion as 
well as within a RPA treatment of the bare propagator 
of collective spin-1 excitations weakly coupled to BCS 
quasiparticles in the spirit of Refs. 85-89. In both cases 
we observe a softening of the q-dependent spin gap as a 
result of switching on higher-order <i-wave gap harmon- 
ics which is rooted in the same mechanism as for the 
single-band Hubbard model. 



V. CONCLUSIONS 

In conclusion, we have examined the effects of higher 
c?-wave gap harmonics induced by underdoping on the dy- 
namical magnetic susceptibility of high-T c cupratcs based 
on a fermiology approach. The calculations are carried 
out for a single-band Hubbard model with an on-site re- 
pulsion treated within the random phase approximation. 
The input parameters for the BCS dispersions are taken 
directly from angle resolved photoemission measurements 
on Bi2212. We find that the inclusion of higher harmon- 



ics decreases the q-dependent spin gap to a degree consis- 
tent with experiments performed on YBCO. This effect is 
robust in that it does not depend on the detailed nature 
of the fermiology model. The downward dispersion of the 
incommensurate peaks is reproduced and shown to move 
down in energy with underdoping. We find a crossover 
from parallel to diagonal incommensuration above the 
resonance energy. However, this effect depends sensi- 
tively on the shape of the BCS dispersions and on the 
details of the fermiology model, here on the assumption 
of an on-site Hubbard residual interaction. With the ad- 
vent of large enough Bi2212 samples for neutron studies 
it will be possible to compare these predictions to exper- 
iments. 
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APPENDIX A: DEGENERACIES OF min k E 2 (q, k) 

The 2-particle dispersion E 2 (q, k) is 27r-periodic in all 
four variables. In the following, when we refer to the 
first Brillouin zone, we have in mind the k-vectors that 
belong to the primitive cell ]— tt, tt] x ]— tt, tt]. For any 
wave vector q, i? 2 (q, k) is invariant under the involutivc 
symmetry transformation 

(k x , k y ) ► (—k x — q x , —k y — q y ), (Ala) 

which can be viewed as a rotation of angle tt about the 
point 

(fcxA) = -(f,f)- (Alb) 

For special values of q the 2-particle dispersion possesses 
additional symmetries. 

At the wave vector q = (tt, tt), E% is left invariant by 
the point-group transformations of the two-dimensional 
square lattice. Together with transformation (Ala) the 
eight elements of the point group form a group of 16 
elements. Hence, the two-dimensional primitive cell in 
k space can be divided into 16 triangular subcells which 
transform into each other under the action of the symme- 
try group. As the 2-particle energy is a continuous and 
bounded function in k, Ei(tt, tt, k) must have at least one 
minimum and one maximum in each subccll. Further- 
more, since the sides of the subcells are mirror axes, the 
gradient of the 2-particle dispersion at q = (tt, tt) van- 
ishes on every corner of the triangles. The degeneracy 
of the (local) minima of Ei(tt, tt, k) in the first Brillouin 
zone is a direct consequence of the symmetry transforma- 
tions. Depending on whether the minima lie in the in- 
terior of the triangles, on the reflection axes, on corners 
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FIG. 7: Panel (a) displays the k dependence of _E 2 (q, k) at 
the q- vector (ji, ir) in the first Brillouin zone ]— n, n] x ]— n, n] 
for the BCS dispersion of Bi2212 at optimal doping (see ta- 
ble I). Panels (b), (c), and (d) are the same as panel (a) 
except for the momentum transfer q = (1.12257T, 1.12257r), 
q = (1.267r,7r), and q = (1.0657T, 0.8057r), respectively. The 
elevation of the contour lines increases with color, red being 
the lowest and violet the highest. The axes of reflection are 
drawn in black. There are one 8-fold degenerate minimum in 
panel (a) , two 2-fold degenerate minima and one 4-fold degen- 
erate minimum in panel (b) , two 4-fold degenerate minima in 
panel (c), and four 2- fold degenerate minima in panel (d). 



which are shared by four triangles or on corners which 
are shared by eight triangles, the (local) minima of the 
2-particle energy are either 16-fold, 8-fold, 4-fold, or 2- 
fold degenerate, respectively. For the BCS dispersion of 
Bi2212 at optimal doping (see table I) the symmetries 
and the degeneracy of the global minimum of -£2(1", tt 5 k) 
are depicted in Fig. 7a. 

When q is chosen to be on either of the two diagonals 



passing through (ir, it), the 2-particle energy is symmetric 
under the reflection about the corresponding diagonal in 
the Brillouin zone. Together with transformation (Ala) 
this yields a symmetry group with four elements. The 
first Brillouin zone decays into four polygonal subcells 
which are related to each other by the action of the sym- 
metry group. Hence, the degeneracy of the (local) min- 
ima is 4-fold, when they arc in the interior of the poly- 
gons, 2-fold, when they lie on a reflection axis and 1-fold, 
when they are on a corner of the subcells. This is illus- 
trated in Fig. 7b for the case when q is on the diagonal 
connecting (— 7r, 7r) and (7r,7r). 

If q is located on the horizontal (vertical) axis passing 
through the commensurate point (71% 7r), E2 is left invari- 
ant under the reflection about the horizontal (vertical) 
axis which goes through the origin of the Brillouin zone. 
The (local) minima arc doubly degenerate, when they lie 
on the reflection axis or on the point (Alb); otherwise 
they are 4-fold degenerate (see Fig. 7c). 

Finally, when q is chosen not to be on any of the above 
symmetry axes, the symmetries of the 2-particle energy 
are reduced to rotation (Ala). The degeneracy of the 
(local) minima is 1-fold, when they lie on the symmetry 
point (Alb), and 2-fold otherwise (see Fig. 7d). 

From the continuity of ^(q, k), now considered as 
a function of q, it follows that the number of local 
minima of the 2-particlc energy in the first Brillouin 
zone for the k's is constant on an open neighborhood 
of q = (ir, 7t). For the BCS dispersions of table I this 
neighborhood includes the area around the antiferromag- 
netic wave vector where most of the spectral weight of 
X^ PA (w,q) is located, as is illustrated in Fig. 7. With 
these hopping and gap parameters the minima of the 
2-particle energy at the commensurate point (it, tt) are 8- 
fold degenerate, since they lie on the magnetic Brillouin 
zone boundary. There is one doubly and one quadru- 
ple degenerate local minimum, and a doubly degenerate 
global minimum, when the momentum transfer lies on a 
diagonal. If q is chosen to be on a vertical or horizontal 
axis passing through the commensurate point (n, 7r), the 
two (local) minima are 4-fold degenerate. Finally, when 
the wave vector q is not on any of the symmetry axes, 
there are four doubly degenerate (local) minima. 
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